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ABSTRACT 

We present high-resolution mid-infrared (MIR) imaging, nuclear spectral energy distributions 
(SEDs) and archival Spitzcr spectra for 22 low-luminosity active galactic nuclei (LLAGN; Lboi ^ 
lQi2 gj,g g-i). Infrared (IR) observations may advance our understanding of the accretion flows in 
LLAGN, the fate of the obscuring torus at low accretion rates, and, perhaps, the star formation 
histories of these objects. However, while comprehensively studied in higher-luminosity Seyferts and 
quasars, the nuclear IR properties of LLAGN have not yet been well-determined. We separate the 
present LLAGN sample into three categories depending on their Eddington ratio and radio emission, 
finding different IR characteristics for each class. (I) At the low-luminosity, low-Eddington ratio (log 
Lboi/LEdd < —4.6) end of the sample, we identify "host-dominated" galaxies with strong polycyclic 
aromatic hydrocarbon bands that may indicate active (circum-)nuclear star formation. (II) Some 
very radio-loud objects are also present at these low Eddington ratios. The IR emission in these 
nuclei is dominated by synchrotron radiation, and some are likely to be unobscured type 2 AGN that 
genuinely lack a broad line region. (HI) At higher Eddington ratios, strong, compact nuclear sources 
are visible in the MIR images. The nuclear SEDs of these galaxies are diverse; some resemble typical 
Seyfert nuclei, while others lack a well-defined MIR "dust bump" . Strong silicate emission is present 
in many of these objects. We speculate that this, together with high ratios of silicate strength to 
hydrogen column density, could suggest optically thin dust and low dust-to-gas ratios, in accordance 
with model predictions that LLAGN do not host a Seyfert-like obscuring torus. We anticipate that 
detailed modelling of the new data and SEDs in terms of accretion disk, jet, radiatively inefficient 
accretion fiow and torus components will provide further insights into the nuclear structures and 
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processes of LLAGN. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - galaxies: photometry — 
infrared: galaxies 



1. INTRODUCTION 

An active galaxy spends only a small fraction of its 
lifetime in a spectacular, highly luminous, qua sar-like 
phase (e.g. [Hopkins & Hernquist 2006; Shanka r" et al.l 
I2009f ). Much more time is spent in a near-quiescent, 
weakly-accreting state, and indeed low-luminosity A GN 
(LLA GN) are found in nearly half of nearby galaxies (jHol 
I2008f ). The properties of the accretion flow and associ- 
ated structures in LLAGN are thought to be quite dif- 
ferent from those in their higher-luminosity counterparts. 
A full description of AGN accretion and lifecycles, then, 
will require understanding the low-luminosity members 
of the class. 

In this paper we investigate the nuclear infrared (IR) 
properties of LLAG^Q. The IR continuum emission from 
an AGN may contain contributions from a variety of pro- 
cesses: thermal emission from a dusty torus surround- 
ing the central supermassive black hole (SMBH), syn- 
chrotron emission from a jet, and - in LLAGN - ther- 
mal emission from a truncated accretion disk (see be- 
low). Dust in a narrow-line region (NLR) or associated 
with the surrounding stellar population may also con- 
tribute. Silicate emission or absorption features at 10 
and 18 /xm contain information about dust geometry and 
heating, and the suite of IR polycyclic aromatic hydro- 
carbon (PAH) bands may trace star formation in the 
nuclear environment. The IR wavelength regime there- 
fore contains much information that could potentially 
advance our understanding of the nature and lifecycle 
of LLAGN. However, IR observations of LLAGN remain 
fairly scarce, especially at the angular resolution neces- 
sary to isolate the weak nuclear emission from that of the 
host galaxy. 

Studies of the spectral properties, luminosities, fuel 
supply etc. of LLAGN indicate that the SMBHs in 
these objects are accretin g at low rates and with low 
radiat ive eflnciencv (^^1^1 120071: IHoI 120091: iTrump et ahl 
120111: lYu et al.l |201l|). These studies suggest that the 
standard optically thick and geometrically thin accre- 
tion disk is truncated in LLAGN, and replaced by 
a geometrically thick, radiatively inefficient accretion 
flow (RIAF: e.g. IRees et all [19821: iN aravan fc Yil [1991 
iNaravan fc McClintockll2008[ ) interior to the disk trun- 
cation radius. In addition, jet emission becomes in- 
creasingly importa nt as the accretion rate decreases 
([Yuan fc Cuill2005| ). The truncation of the thin disk 
should shift its thermal emission peak towards longer 
wavelengths, leading to the suggestion that the mid-IR 
(MIR) peaks or excesses observed in the SEDs of some 
LLAGN are the equi valent of the "big blue bump" of 
Seyf erts and quasars (|Holl2008l ). Indeed. [Nemmen et al.l 
(pOlli) present models in which the bulk of the nuclear 

^ For the purposes of this paper, we loosely define a "low- 
luminosity" AGN as a low-ionisation nuclear emission region 
(LINER) or Seyfert galaxy with Lbol below about 10*^ erg s~-^. 
For comparison, Sgr A* has Lboi 10^'' erg s~^, the dwarf Seyfert 
NGC 4395 has Lboi ~ 10'"' erg and the luminous quasar 

PDS 456, Lbol ~ 10*^ erg s"!. 



luminosity at ^1-10 fim can come from a truncated thin 
disk. However, the authors note that the available, large- 
aperture IR data are not sufficient to constrain the pres- 
ence or properties of a truncated disk in the objects mod- 
eled. IR photometry, at higher angular resolution than 
is generally available in the literature, may therefore aid 
our understanding of the accretion physics of LLAGN. 

In higher-luminosity AGN, the difference between 
types 1 and 2 is explained at least to first order by the 
presence of a toroidal cloud of dust and gas obscuring the 
AGN from certain viewin g directions whil e permitting a 
direct view from others (IAntonuccilll993D . An obvious 
question, then, is how (if at all) do LLAGN fit into this 
unified framework? 

Some models of the torus explain its existence through 
infiows of gas from larger scales, often inv oking the ef- 
fects of nucle ar star clusters or disks (e.g . 'Wada et al.l 
2009; Schartm a^ et al.ll2010l:lHopkins et al. 2011). Dur- 
ing periods when little material is reaching the centre of 
the galaxy, the to rus may become thin and transpar- 
ent (jVollmer et al.l 12008). Conve rsely, the torus inay be 
part of a dusty, outfi owing wind ("Konigl fc Kartid ll994l : 
lElitzur fc Shlosman 2006: Dorodnitsvn et al. 2011). The 
disk wind model of lElitzur fc Shlosmanl (2006) predicts 
that below Lboi ~ erg s~^, accretion onto the 

black hole can no longer sustain the outflow necessary 
to obscure the nucleus. In either case, low-luminosity 
AGN may s how little nuclear obscurati on and dust emis- 
sion. In the lElitzur fc Shlosmanl (|2006D model, very low- 
luminosity AG N are also expected to lack a broad line 
regio n (see also lNicastrol[2000t lLaoril2003t lEhtzur fc Hoi 
[20091 ). 

Observationally, there are indications that LLAGN do 
tend to have unobscured nuclei. For instance, several 
authors have detected nuclear UV/o ptical point sources 
in LINERs of both types 1 and 2 (iPogge et al.l 120001: 
[Chiaberge et al.l l2005t iMaoz et al.l 120051 ) . In a study 
of 86 nearby galaxi es with point l ike X -ray nuclei and 
Lx < 10''2 erg s'^, iZhang et al.l ([2009!) find iVn to be 
correl ated with X-ray luminosity. iBrightman fc Nandral 
(|2011f ) also find a decline in obscuration at low lumi- 
nosities. Detailed X-ray observations of some LLAGN 
have shown that the Fe Ka line, a signature of X-ray 
reprocessing by cool material, is w eak in these nuclei 
(jPtak et al.ll2004l: [Binder et al.i[2009t ). 

However, this is not the case for all objects. The 
detect ion of broad Ha i n polarised light in 3/14 LIN- 
ERs (iBarth et al.l (1999) suggests that these objects 
host dust-obscured AGN, as do the biconic al ionization 
cones observed in a handful of LINERs (jPogge et al.l 
|2000( ). X-ray studies of single objects show that some 
LLAGN do have substantial a bsorbing columns (e.g. 
NGC 4261; IZezas et al.l [200l, and the fraction of 
such o bje cts may be s ignificant (jConzalez-Martm et al.l 
l2009al ). iSturm et al.l (2006) present average spectra 
of type 1 and 2 LINERs which suggest an extra hot 
dust component present in the type Is relative to the 
type 2s, similar to the SED s of "conv entional" Seyferts 
((Alonso-Hcrrero et al.,.2003:,Deo et al.i,2009;.Prieto et all 
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[Mot iRamos Almeida et al.|[20Tll ). The role of dust in 
LLAGN, and in particular whether there exists a lumi- 
nosity, accretion rate, or other property at which the 
torus ceases to exist, remains unclear. A search for the 
IR signatures of the torus - thermal emission and silicate 
emission and absorption features - promises a better un- 
derstanding of these issues. 

Finally, IR observations may give information about 
the star formation histories of LLAGN. Based on op- 
tical spectroscopy, LINE R stellar p opulations ar e gen- 
erally thought to be ol d (|Ho et al.i r2003; Kcwlc v et al.l 
I2006t [Zhang et al]|2008D . although intermediate- age pop- 
ulations have been det e cted in a minority of cases 
([Gonzalez Delgado et al.l I2004D . However, the opti- 
cal regime is not ideally suited for discovering very 
young stars few Myr) or star formation activ- 
ity. Small-aperture UV spectroscopy of stellar absorp- 
tion lines can be a useful tec hnique (|Maoz et al.l 119981 : 
iStorchi-Bergmann et al.|[2005[ ). but is vulnerable to ex- 
tinction. Althoug h not always an unambig uous st ar for- 
matio n indicator ([Peeters et al.|[200^ iDia z-Santo s et all 
I2010f) . observations of PAH emission in LLAGN could 
prove to be a complementary method of tracing young 
stellar populations and/or ongoing star formation (e.g. 
[Shapiro et al.l [20100 . IR spectroscopy can also probe 
older stellar populations, which may be important in ex- 
pla ining the energy budgets of LLAGN fEra cleous et al.l 
I2010a|) . by locating spatially exten ded silicate emission 
features from AG B star envelopes ([Bressan et al.l [20061 : 
iBuson et al.l[200l . 

At the time of writing, the IR emission of LLAGN 
remains comparatively unexplored, particularly on spa- 
tial scales < 1". For many years following early , 
Jower-rcsolution MIR work fe.g. [ Ro che fc Aitken lll985l: 
Lawrence jrt_aL 198 5; Will ner et al, Tri985: Cizdzi el et al.l 
1981 iRoche etm fi99it IKnapp et al.l [T9p), high- 
resolution imaging at A ~ 10 fim existed only 
for a handful of r e latively bright , well - known objects 



jCharv et al.' '2000^ 'Gr ossan et all [20011: Ferl man et al.' 
2001; Whvsong & An tonuccil l2004t [Mason e t al. 2007; 
Horst et al.H2008l : iRadomski et a l. 2008; Rcunanc n et al.l 
More recentlv. I Asmus et~a\. (2011.) have presented 



ground-based MIR imaging of a num ber of LLAGN, with 
seven new detections. In addition, Ivan der Wolk et all 
([2010[ ) detected 10 radio galaxies in their high-resolution 
MIR images, some of which can be considered LLAGN. 
Some subarcsecond-resolution NIR imaging is also avail- 
able in the literature, but it has rarely been considered in 
the context of the multi- wavelength emission of LLAGN. 
Published LLAGN SED compilations contain little or no 
high-resolution IR data ('Ho"1999l :Ftak et al.|[200llMaozl 
[2007; Eraclcous ct al. 2010b). 

With the aim of illustrating the overall nuclear IR 
properties of LLAGN, we have acquired new, ground- 
based MIR "snapshot" imaging of 20 IR-faint nuclei. 
Acquired with Michelle and T-ReCS on the Gemini tele- 
scopes, the angular resolution of the observations is ap- 
proximately 0".35 at A ~10 /im. We have also obtained 
seeing-limited images at 3 - 5 /im of five objects. To 
the new data, we add published, high-resolution IR pho- 
tometry of these and a further two nuclei. We combine 
these data with published measurements at other wave- 
lengths to produce nuclear spectral energy distributions 
(SEDs) for these objects. In terms of high spatial reso- 



lution IR data, these SEDs are by far the most detailed 
yet available. Finally, we also present archival Spitzer 
low-resolution spectroscopy for the 18/22 galaxies with 
available data. 

In [21 we discuss the galaxy sample studied in this pa- 
per. SJ3] describes the reduction of the MIR and NIR ob- 
servations and the resulting photometric measurements. 
In ij4.1l we present the images and use the MIR/X-ray re- 
lation to investigate the nature of the structures revealed 
in them. i i4.2l shows the IRS spectra of the LLAGN, and 
the full, high-resolution radio - X-ray SEDs are presented 
in ij4.3l (with a detailed view of the 1-20 /im region in 
Appendix [^. The processes that may give rise to the IR 
emission in the different types of LLAGN are discussed 
in fj5] We anticipate that the new IR data and SEDs 
will provide valuable constraints on models of the nu- 
clear processes in LLAGN, but we defer such detailed 
modelling to future work. 

2. THE LLAGN SAMPLE 

The 22 LLAGN studied in this paper have L2 -10 k cV < 
3 X lO''^ erg (corrected for extinction; see i )4.3p and 
D < 36 M pc (1"< 170 pc). Assuming Lboi/Lx ^ 16 
([Holl2009H^ . this X-ray luminosity corresponds to Lboi ~ 
5 X 10"^^ erg s~^, roughl y the threshold below which 
lElitzur fc Shlosmanl ([20061 ) predict the disappearance of 
the t orus in LLAGN. Spectral classifications were taken 
from [Ho et "all ([1997aD . Most of the objects are LINERs, 
but the sample also contains a few low-luminosity Seyfert 
galaxies. Basic information for all of the sources is given 
in Table [H 

All of the galaxies we observed had old, low- 
resolution IR photometry available i n the literature (e.g. 
[Willner et allll985l: iDevereux et~all I1987D . and for this 
pilot program we favored objects which had been de- 
tected in the past. The sample as a whole may there- 
fore be biased towards galaxies which are bright in the 
IR. The LINERs in this sample are, however, "IR-faint" 
in terms of their FIR/optical luminosity ratios. Values 
of JjFiji/hB range from 0.1-5 for this sample, com- 
pared with median values of I^fir/^b of 106 and 1.8 
for th e IR-bright and IR-faint objects in iSturm et al.l 
([2006( 1. respectively. IR-bright and IR-faint LINERs 
are cleanly separated in several MIR diagnostic dia- 
grams, and inay re present entirely different phenomena 
([Sturm et al] 120061) . Although many IR-bright LINERs 
show the compac t hard X-ray sour ces characteristic of 
AGN activity fe.g. lDudik etalEoO^ . the optical LINER 
emission in those objects may not be directly related to 
the accretion onto the central black hole (Ho 2008) . 

There is evidence that each of the galaxies in the 
sample genuinely hosts an actively accreting, supermas- 
sive black hole. This evidence takes the form of an 
unresolved hard X-ray point source, optical/UV vari- 
ability, high-ionisation MIR emission lines, and other 
phenomena. The specific evidence for each object is 
given in Table [TJ We estimate the Eddington ratios, 
Lboi/LEdd, of the AG Ns using the stellar velocity dis- 
persions of iHo et al.[ (p309), the Mbh - ^ relation of 

* All bolometric luminosities (and, therefore, Eddington ratios) 
given in this paper are derived using this conversion factor unless 
otherwise stated. While a convenient tool, X-ray-based bolometric 
correc tions are still a matter of some debate; see, for example, 
Eracleous et al.1 l|2010bl l. 
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Figure 1. Absorption-corrected hard X-ray luminosity ys 
Eddington ratio for the LLAGN sample. See !j2]for definitions 
of the three LLAGN classes. 



iTremaine et al.l (|2002f ) and X-ray based bolometric lu- 
minosities. As expected, the Eddington ratios are l ow: 
6 X 10 '''' < Lboi/LEdd < 4 X 10"'*. For comparison, IHoI 
(|2009D finds a median Lboi/LEdd ~ 10"^ for the LINERs 
and Lboi/LEdd ~ lO""* for the Seyferts in the Palomar 
sample. 

In 12] onwards we divide the LLAGN into three sets for 
the purpose of discussion: 

L At log Lboi/LEdd < —4.6, with PAH-dominated spec- 
tra and lacking a strong nuclear point source, are the 
"host-dominated, low-Eddington ratio" sources. 

II. Also at log Lboi/LEdd < —4.6 we highlight a group 
of LLAGN with s trong radio emission (log v L,j{'a GHz) 
/ Lx > -2.7; cf iTerashima fc WilsonI |20Q1 and well- 
sampled SEDs that allow strong constraints to be placed 
on the origin of their IR emission. These are the "radio- 
loud, low-Eddington ratio" objects. 

III. Galaxies with log Lboi/LEdd > —4.6 (the Edding- 
ton ratio of NGC 3718, which has the lowest value of 
the objects - aside from NGG 4486/M87 - with a com- 
pact nucleus (|4.ip ). form the "high Eddington- ratio" 
category. 

We choose to discuss the LLAGN in terms of Edding- 
ton ratio as opposed to some other quantity (e.g. lumi- 
nosity, morphology, spectral features) as there are indi- 
cations that accretion rate is the fundamental parameter 
governiiig the properties of these objects (e.gjy uan 2002 
Holl200llElitzur fc Holl2009tlTrumD et al.il2011tlYu et al] 
201 It ) . As shown in Figure [1] Eddington ratio is loosely 
correlated with hard X-ray luminosity in this sample. 
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Optical spectroscopic classification as given by IHo et aLl I ll997al ). except for NGC 1097. IPhillips et al.l 119841) find 
that NGC 109 7 has a "c lassic" LINER spectrum, and the detection of broad components to both the Ha and H/3 lines 
fStorchi-B crgm ann et aI] [T993^ suggest that it should be classified as a type 1. 

^ X = Lbol/LEdd 

" Rx = log uL„(5 GHz) / Lx 

X-ray: Hard X-ray point source(s). Ha: broad Ha line. Pol Ha: broad Ha detected in polarized light. Radio: high 
brightness-temperature radio nucleus and/or radio jet. UV: UV variability. [Ne V]: [NeV] 24 fim line. NLR: ionization 
cone (s). Th is is not intended to be an exhaustive list for each object. 

1: IHo et al. (1997b). 2: Earth ct al. (1999). 3: ISatvapal et al.l 12001'). 4: ^Gonzalez-Martin et al.l SWO^). 5: IDudik et al.l 
(120091). 6: Poggc ct al. (2000). 7: Satvapal et al. (2005). 8: Fa lcke et al.l (120001) . 9: ITera shima fc Wilson (2003). 10: 
Nagar et al (2005). 11: Maoz et _al. (2005); 12: Flohic ct al. (20 0^). 13: IHo et alj (|2001l). 1 4: Laine & Beck (2008). 15: 
Dudik et all (|2007|). 16: iNemmen et al.l (J200Q). 17: Storchi-Bc rgmann et al.l ([1993 ). 18: IGonzalez-Martfn et al.i (l2009h). 
' I: Host- dominated, low-Eddington ratio galaxy (log L^oi/LEdd < —4.6). II: Radio-loud, low-Eddington ratio galaxyflog 
Lbol/LEdd < -4.6, log uL^{5 GHz) / Lx > -2.7). Ill: High-Eddington ratio galaxy (log Lboi/Lfidd > -4.6). See ^lov 
details. 
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3. OBSERVATIONS AND DATA REDUCTION 
3.1. Gemini Mid-Infrared Imaging 

MIR imaging; observatio ns were performed using T- 
ReCS (jTelesco et a l.l ll998D on the 8.1 m Gemini South 
telescope and Michelle (jGlasse et al.lll997D on the Gem- 
ini North telescope. Both T-ReCS and Michelle use a 
Raytheon 320 x 240 pixel Si:As IBC array. This provides 
a plate scale of 0".09 and 0".10 per pixel respectively, 
corresponding to a field of view (FOV) of 28" x 21" 
for T-ReCS and 32" x 24" for Michelle. The standard 
chop-nod technique was used to remove time- variable sky 
background, telescope emission and 1/f detector noise. 
In all observations the chop throw was 15", the chop 
angle was 30 ° E of N and the telescope was nodded ap- 
proximately every 40 s. 

The Si-2 (8.7 ^im, AA = 0.8 ^im, 50% cut-on/off) filter 
in T-ReCS provides the best combination of sensitivity 
and angular resolution, so this filter was used for the 
majority of the T-ReCS observations. NGC 1052 was 
also observed with T-ReCS' Qa (18.3 /zm, AA = 1.5 /um, 
50% cut-on/off) filter. Most of the Michelle observations 
were taken using the sensitive semi-broad N' filter (11.2 
fim, AA = 2.4 ^m, 50% cut-on/off). Some of the brighter 
galaxies were also observed with Michelle in one or more 
of the Si-2 (8.8 /xm, AA = 0.9 ^m, 50% cut-on/off), Si- 
4 (10.3 ^m, AA = 1.0 fim, 50% cut-on/off), Si-6 (12.5 
^m, AA = 1.2 ^m, 50% cut-on/off) and Qa (18.1 ^m, 
AA = 1.9 fim, 50% cut-on/off) filters. For galaxies with 
existing, published high-resolution MIR photometry, the 
Michelle filters were chosen to complement the available 
data. A summary of the observations, which were taken 
in queue mod^EI, is given in Table [51 

The data were reduced using the Gemini IRAF pack- 
age and our own IDL routines. The difference for each 
chopped pair was calculated and the nod-sets then dif- 
ferenced and combined to create a single image. During 
this process all nods were examined for high background 
that could indicate the presence of clouds or high water 
vapor, but no data needed to be removed for this reason. 
In a few cases where AB pairs showed unusually high 
electronic noise, those pairs were reje cted. 

Observations of Gohen standards (jGohen et al.lll999h 
were obtained immediately before and/or after each 
galaxy observation using the same instrument configu- 
ration. A Moffat function with two free parameters, 
FWHM and best described the delivered PSF (cf 
iRadomski et al.ll2008[ ). The galaxy images were fiux cal- 
ibrated using standard star brightnesses appropriate for 
the instrument, filter, airmass and standard star SEEQ. 
At this stage an artificial source of known flux (10 mJy) 
was introduced into the calibrated images, with the same 
FWHM and /? parameter as the associated standard star. 
This source was used to quantify changes to the signal-to- 
noise ratio (S/N), galaxy brightness and FWHM during 
the subsequent reduction steps, and to provide an esti- 
mate of the error on the photometric measurements. 

Vertical and/or horizontal electronic patterns in the 
images were then removed, if necessary. This was 

8 Program IDs; GN-2007A-Q-93, GS-2007B-Q-203, GN-2008A- 
Q-43, GN-2011A-Q-55 

Calculat ed using the web form available at 
|http: / / www .gemini.edu/?q=node/11246| 



achieved using the "miclean" task in the Gemini IRAF 
package, which sequentially models and removes elec- 
tronic noise and background structures. Finally, in the 
cases where the galaxy was observed more than once, the 
data were combined. 

Two different photometric measurements were made 
for each galaxy (Table [3]). First, the flux in a circu- 
lar aperture of 2" diameter was measured. Second, the 
standard star was scaled to the centroid of the galaxy 
emission and then the flux measured in a 2" aperture 
around the scaled star. We refer to these methods as 
"aperture" and "PSF-scaling" photometry, respectively, 
in the remainder of the paper. The PSF-scaling photom- 
etry represents the maximum likely contribution from 
an unresolved nuclear source, while the difference be- 
tween the aperture and PSF-scaling measurements gives 
the approximate contribution of extended emission sur- 
rounding the central source. Various authors have fur- 
ther refined the PSF-scaling technique by adjusting the 
scaling of the standard star to account for the underly- 
ing host gal axy emission that undoubtedly contributes at 
the nucleus fRadomski et al. 2003|; iLevenson et al.ll2009l : 
Ramos Almeida et al. 2009, 201 1)|). However, because of 
the low S/N in some of the images, and to provide a 
uniform set of measurements, we did not attempt this. 
The above groups typically find that in Seyfert galaxies, 
70 - 100% of the simple PSF-scaling flux in the N band 
comes from the pointlike nucleus. We therefore expect 
the PSF-scaling photometry to provide a good estimate 
of the AGN brightness in LLAGN dominated by nuclear 
point sources. This will not be true for objects with very 
faint nuclei embedded in copious extended emission. 

The uncertainty in the photometry was estimated in 
the following manner. First, the time-variability of the 
sky transparency was estimated from the variation in 
the signal from the standard stars, found to be ~6% at 
A 10 /im and ~10% at A ~ 18 /zm. Second, the dif- 
ference between the initial and final fluxes of the artifi- 
cial source was calculated for each galaxy. Third, for the 
PSF-scaling photometry, an estimate of the error induced 
by a variable PSF was obtained by cross-calibrating the 
standard stars observed on a single night. This error was 
found to be ~ 13%. The total uncertainly was calculated 
by adding in quadrature these individual contributions. 
The pixel-pixel variations in the residual sky background 
contribute a negligible amount to the photometric errors. 

All galaxies but one, NGC 4374, were detected in the 
MIR imaging. We use the observations of NGC 4438, a 
faint object observed on the same night, to estimate a 
3(7 upper limit on the 8.7 /xm flux of NGC 4374. To 
complete the set of consistently obtained photometric 
measurements, we al so re- analysed archi val MIR imaging 
data for NGC 1097 (iMason et all [2001 and NGC 4486 
(jPerlmanet al.l[200l . 

3.2. Gemini Near- Infrared Imaging 

NIR imaging observations (GN-2011A-Q-55) were per- 
formed in queue mode on 20110204 using th e Near In- 
fraRe d Imager and Spectrometer (NIRI; Ho dapp et al.l 
I2003D on Gemini North. The NIR observations were ob- 
tained for a subset of the brighter galaxies in the sample 
and were intended to fill in parts of the SEDs lacking in 
published data. 

NIRI uses a 1024 x 1024 ALADDIN InSb array. With 
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Table 2 

Summary of new MIR observations 



Galaxy 


Instrument 


Date 


Filter 


Integration Time^ 
(s) 


NGC 1052 


T-ReCS 


20070817 


Si2 


927 




T-ReCS 


20070820 


Qa 


927 


NGC 3031 


Miehelle 


20110203 


Si2 


311 




Michelle 


20110203 


Si4 


309 




Michelle 


20110203 


Si6 


329 




Michelle 


20110204 


Qa 


324 


NGC 3166 


T-ReCS 


20071202 


Si2 


927 


NGC 3169 


T-ReCS 


20071203, 20071206 


Si2 


1853 


NGC 3718 


Michelle 


20081107 


N' 


855 


NGC 3998 


Michelle 


20110123 


Si2 


311 




Michelle 


20110123 


Si4 


309 




Michelle 


20081206 


Si6 


753 




Michelle 


20110123 


Qa 


324 


NGC 4111 


Michelle 


20081206 


N' 


743 


NGC 4258 


Michelle 


20110203 


Si2 


311 




Michelle 


20110203 


Si4 


309 




Michelle 


20110203 


Si6 


329 


INGC 42b L 


T-ReCS 


20080103 


Si2 


927 


NGC 4278 


Michelle 


20080616 


N 


525 


JNGC 4374 


T-ReCS 


20080220 


Si2 


927 


NGC 4438 


T-ReCS 


20080220, 20080704 


Si2 


1853 


NGC 4457 


T-ReCS 


20080327 


Si2 


927 


NGC 4579 


Michelle 


20110205 


Si2 


311 




Michelle 


20110123 


Qa 


648 


NGC 4594 


T-RoCS 


20080502 


Si2 


927 


NGC 4736 


Michelle 


20070415 


N' 


612 




Michelle 


20070416 


Qa 


512 


NGC 5005 


Michelle 


20081206 


N' 


900 


NGC 5033 


Michelle 


20080616 


N' 


941 


NGC 5363 


Michelle 


20081107 


N' 


900 


NGC 7479 


Michelle 


20080701 


Si2 


933 



On-source integration time, excluding any rejected nods 

estimate a Scr upper limit on the M' flux density of 
NGC 4486. 

3.3. Spitzer IRS spectra 

Archival low-resolution spectra taken with the 
Spitzer Space Tel escope's InfraRed Spectrograph (IRS; 
iHouck et al]|2004D are available for 18 of the 22 LLAGN 
and were downloaded from the Spitzer science archive^. 
Starting from the pipeline-processed "bed" data, the files 
were background-subtracted either by subtracting spec- 
tra in different nod positions, or by subtracting sky ob- 
served in one order when the object was being observed in 
the other order, depending on the available observations. 
One-dimensional spectra were extracted using the point 
source extraction algorithm provided by the SPICE soft- 
ware package, which provides a wavelength-dependent 
aperture of 7.2" at 6 fim. A few objects were observed 
in spectral mapping mode; for these only the spectrum 
at the nuclear position was extracted. Inspection of the 
spatial profiles of the spectra showed that many of the 
galaxies are extended at Spitzer's angular resolution, to 
varying degrees. The slit loss correction applied by the 
SPICE point source flux calibration is appropriate for 
an unresolved source and cannot account for wavelength- 
dependent extended structure. This may affect the accu- 
racy of the flux calibration and overall continuum slope 
in the spectra of the extended objects in this sample. 
iGaUimore eFal] (|2010D find errors of ~10% from this ef- 
fect in their analysis of the 12 AGN sample (contain- 
ing HII region galaxies, LINERs and Seyferts). Nonethe- 



the f/32 camera, used to avoid saturation on the bright 
sky background, this provides a plate scale of 0.022" per 
pixel and FOV of 22.5" x 22.5". Observations were taken 
through the L' (3.9 /zm, AA = 0.4 /xm, 50% cut-on/off 
) and M' (4.8 ^m, AA = 0.1 fim, 50% cut-on/off ) fil- 
ters. Total integration times of 180 sec (L') and 225 sec 
(M') were used. The galaxy images were acquired in an 
ABB'A' pattern, nodding approximately 50" to blank 
sky. 

The data were reduced using the Gemini IRAF package 
and standard NIR imaging procedures. The images were 
corrected for nonlinearity, sky-subtracted, and divided 
by a flat fleld constructed from sky frames. A bad pixel 
mask was also applied. The dithered observations were 
then shifted to a common position and averaged, and flux 
calibration was achieved using standard stars observed 
imme diately before or after each galaxy (jLeggett et al.l 
I2OOI . 

As with the MIR observations (> j3.1|) . an artiflcial 
source was introduced at this stage and both aperture 
and PSF-scaling photometry performed (Table 3]). For 
NGC 3031 and NGC 4258 at M', the FWHM of the stan- 
dard star is larger than that of the AGN, causing the 
PSF-scaling flux to be higher than that given by the 
aperture photometry. Therefore no PSF-scaling fiux is 
quoted for these objects. The uncertainty in the pho- 
tometry was calculated in the same way as for the MIR 
photometry ( ij3.ip . 

NGC 4486 was not detected in these short M' band 
exposures. We use the observations of NGC 4258 to 
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Table 3 

MIR photometric results 



Galaxy 


Filter 




Fpsf'^ 


FWHMgai'^ 


FWHMstd'* 






(mJy) 


(mJy) 


// 


// 


JNGC lOoi 


Si2 


62 ± 4 


42 it 7 


0.42 


0.39 




Qa 


165 ± 17 




0.75 


0.67* 


NGC 1097^ 


Si5 


42 ± 3 


39 it 6 


0.43 


0.41 




Qa 


64 ± 4 


63 it 9 


0.54 


0.52 


M Q n Q 1 


DIZ 


69 it 5 


67 it 10 


0.37 


0.36 




Dl4 


145 it 10 


140 it 21 


0.39 


0.35 




olD 


1 on -1- Q 
Izy it o 


lie; J- 17 
ilo it 11 


U.4o 


u.oy 




Qa 


zoi ± z6 


1 07 _1_ QO 

197 ± oZ 


0.54 


0.53 




Q:o 

bli 


1 1 

1 ± I 


n 1 -1 
O ± 1 




0.27 


M Q 1 n 




13 it 2 


9 it 2 


0.50 


0.46 




IN 


21 it 2 


14 it 3 


0.35 


0.35 


MO/"' Qfino 

inLjC oyyo 


a;o 

bli 


27 it 2 


26 it 4 


0.42 


0.33 




bl4 


57 it 5 


42 it 8 


0.39 


0.34 




Si6 


70 it 5 


52 it 8 


0.54 


0.46 




Qa 


108 it 11 


96 it 16 


0.59* 


0.53 


NGC 4111 


N' 


10 it 2 


4 it 2 




0.35 


NGC 4258 


Si2 


64 it 4 


64 it 10 


0.32 


0.31 




Si4 


87 it 6 


70 it 11 


0.33 


0.32 




Si6 


114 it 8 


108 it 16 


0.37 


0.37 


NGC 4261 


Si2 


5 it 1 


3 it 1 




0.48 


NGC 4278 


N 


5 it 1 


2 it 1 




0.35 


NGC 4374 


Si2 


< 3 






0.32 


NGC 4438 


Si2 


11 it 1 


7 it 2 


0.33 


0.32 


NGC 4457 


Si2 


12 it 1 


6 it 1 


0.47 


0.27 


NGC 4486S 


N 


17 it 1 


16 it 3 


0.46 


0.46 


NGC 4579 


Si2 


32 it 2 


28 it 4 


0.29 


0.24 




Qa 


112 it 12 


97 it 16 


0.58 


0.52 


NGC 4594 


Si2 


8 it 2 


2 it 1 




0.38 


NGC 4736 


N' 


33 it 4 


11 it 4 


0.56 


0.48 




Qa 


101 it 11 


43 it 8 




0.65 


NGC 5005 


N' 


26 it 2 


4 it 1 


0.65 


0.35 


NGC 5033 


N' 


22 it 2 


15 it 3 


0.35 


0.35 


NGC 5363 


N' 


11 it 3 


3 it 2 


0.71 


0.35 


NGC 7479 


Si2 


158 it 10 


157 it 23 


0.40 


0.40 



Photometry in a 2" diameter aperture 
^ PSF-scaling photometry; see text 
" FWHM of nucleus 
<* FWHM of standard star 

For NGC 1052 only, the standard star image in the Qa filter contains a 
prominent diffraction ring. This means that the Moffat function used to 
model the PSF for the PSF-scaling photometry does not provide a good rep- 
resentation of the flux profile in this case. 

* Gaussian FWHM; Moffat function did not converge for this source 

s Re-analys i s of data originally published by IPerlman et al.l 120011) . 
IMason et aT] 1(20071) 



Table 4 

NIR photometric results 



Galaxy 


Filter 






FWHMgai" 


FWHM. 






(mJy) 


(mJy) 


(") 


(") 


NGC3031 


L' 


74 it 2 


56 it 5 


0.42 


0.40 




M' 


56 it 4 




0.38 


0.40 


NGC3998 


L' 


24 it 2 


11 it 2 


0.48 


0.27 




M' 


18 it 3 


16 it 3 


0.39 


0.38 


NGC4258 


L' 


29 it 2 


25 ± 3 


0.47 


0.42 




M' 


23 it 2 




0.28 


0.54 


NGC4486 


L' 


8 it 1 


8 it 1 


0.43 


0.35 




M' 


< 8 






0.50 


NGC 4579" 


L' 


17 it 2 


13 it 2 


0.65 


0.35 



^ Photometry in a 2" diameter aperture 
*' PSF-scaling photometry; see text 

FWHM of nucleus 

FWHM of standard star 
° The M' images of NGC 4579 showed a strong background gradient 
that made reliable photometry impossible 
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less, we expect the IRS spectra presented here to provide 
a useful overview of the general IR spectral characteris- 
tics of the LLAGN. 

4. RESULTS 

4.1. IR Morphology and the MIR/X-ray Relation 

The MIR images are presented in Figures [2] - [5] For 
comparison, archival HST images in various optical/UV 
filters are also shown, where available. The MIR filters 
used in this study have moderate bandwidths and en- 
compass a number of possible emission and absorption 
features. Perhaps most significantly, the Si-2 and N' fil- 
ters may contain PAH emission, while the Si-4 and Qa 
filters can be affected by silicate emission/absorption, if 
present. The [Ne II] fine structure line lies within the 
bandpass of the Si-6 filter. However, clear differences 
in morphology are not seen in objects observed in more 
than one filter and only one filter per object is shown in 
Figures [2] -in 

The images show a wide range of MIR morphologies. 
Some of the nuclei are Seyfert-like, in the sense that they 
are dominated by compact sources with FWHM compa- 
rable to that of the corresponding standard star. All of 
the galaxies observed at 3 - 5 //m fall into this category 
and they all show strong, compact nuclei at those wave- 
lengths as well. Others have a weaker nuclear MIR source 
embedded in a substantial amount of diffuse, extended 
emission. In some cases, notably NGC 5005 (Figure [2]), 
the MIR emission extends over several arcseconds (>100 
pc). 

In the type 1.5 Seyfert NGC 5033 (Figure [5|), faint 
emission regions are detected to the north and south 
of the nucleus, approximately coincident with the [OIII] 
emissio n observed in this obj ect at lower angular reso- 
lution (|Mediavilla et al.ll2005( ). This is likely one of the 
relatively few Seyferts in which the narrow-line region 
(NLR) has been resolved in the MIR. The 18 /im im- 
age of NGC 1052 (Figure |4]) shows faint emission along 
the radi o jet axis, also de tected in optical Ha-|-[N II] 
images (jPogge et al.l l200"ol ) . However, this image does 
show some low-level trefoil structure and the slight ex- 
tension observed in the 8.7 /im image is offset by about 
10° from the 18 /im extended emission, so detection of 
the NLR/jet in NGC 1052 should be regarded as tenta- 
tive. 

4.1.1. Trends m morphology with Eddmgton 
ratio /luminosity 

Figures [5] - [S] demonstrate differences in morphology 
between objects with low and high Eddington ratios. 
Strong nuclear point sources are common, although not 
universal, in the high-Lboi/LEdd objects. Apart from 
NGC 4486 (which is pointlike) and NGC 4374 (which we 
did not detect), all of the objects in the low-Lboi/LEdd 
categories have nuclei that are much weaker compared to 
the surrounding diffuse emission. As discussed in iJTJ the 
presence or absence of a strong nuclear IR source could 
be related to the presence or absence of a dusty torus. 
However, as illustrated by Figure [l] low-Lboi/LEdd ob- 
jects tend also to have relatively low X-ray luminosities, 

^ A reduced IRS spectrum of NGC 3998, including data from the 
SL an d LH modules, was kindly provided by E. Sturm l|Sturm et al.l 
[2001) . 



which could mean that a nuclear point source is simply 
too faint to detect and/or is overwhelmed by host galaxy 
emission in these objects. We now quantitatively inves- 
tigate this point. 

To roughly quantify the extent to which each galaxy 
is dominated by a point source, we take the ratio of the 
PSF-scaling photometry and the 2" aperture photometry 
(|j3l Table [3|). The resulting values range from 0.15 in 
NGC 5005 to 1.0 in NGC 4258. We refer to galaxies 
with Fpsp/F^ii > 0.65 as "point-source-dominated", or 
"compact" , and the remainder of the objects as "diffuse" . 
This division is arbitrary, but it corresponds quite well 
to our subjective, visual classification of objects having 
a prominent nuclear source. 

The well-known MIR/X-ray relation can be used to 
predict whether the torus should be detectable in this 
shallow imaging survey. In Seyfert galaxies, high-energy 
photons are reprocessed into IR radiation by dust in 
the immediate surroundings of the active nucleus. As 
a result of this interaction, a fairly tight correlation 
has been shown to exist between the MIR and hard 
X-ray emission in these objects (with the X-rays es- 
sentially acting as a proxy fo r optical /UV em i ssion ; 
Krabbc et al. 2001: iLutz et al.1 12 004: Ho rst et alJ 1 20061 



2008.: ,Ramos Almeida et al 



2007t_iGandhi et al.' '200£ 



Hardcastle et al. 2009; L evenson et al . 2009; Honig et alJ 
2010"). As the SEDs of LLAGN may be deficient ^in UV 
and optical radiation (Ho 1999, but see also Eracleous 
et al. 2010), the Seyfert-based MIR/X-ray relation may 
overestimate the reprocessed MIR luminosity expected 
from lower- luminosity objects. However, we use the re- 
lation as a starting point to assess the detectability of 
the torus. We fit the combined samp le of the "well- 
resolved" objects in 'Gandhi et al.' (2009") and the "PSF- 
fitting" photometry of .Levenson et al. (2009[ ). measure- 
ments which are intended to most accurately represent 
the nuclear emission of Seyferts and quasars, finding 
log Lx — 0.88 log Lmir + 4.75. The fit is then extrapo- 
lated to the luminosities of the LLAGN. 

Figure [6] shows that the compact, relatively high- 
luminosity, high-Lboi/IjEdd sources are predominantly 
those in which the torus is expected to be detectable 
in this imaging survejH- That we observe strong, unre- 
solved MIR emission in the objects in which we would 
expect to detect the torus is consistent with most or all 
of that emission arising in the torus. Other AGN-related 
emission mechanisms may well be important, though: 
we also find that NGC 4486 and NGC 3718, in which 
the torus should not be detectable, host strong, point- 
like nuclei. We may expect that the IR emission from 
those point-source-dominated objects is associated with 
the AGN rather than the host galaxy, but contains a sig- 
nificant contribution from sources other than the torus. 
We return to this point in i i4.3.2l 

Conversely, the MIR/X-ray relation predicts that torus 
emission in most of the diffuse objects should be too 
faint to be detected. While the objects with diffuse 
morphologies mainly come from the low-Eddington ra- 
tio categories, the lack of a dominant pointlike nucleus 

This simple calculation does not take into account the diffi- 
culty of isolating a faint torus surrounded by relatively bright host 
galaxy emission, but for the purposes of interpreting Figure|6l that 
is only important for galaxies just above the threshold. 
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NGC41 1 1 (N') 


NGC41 1 1 (F457M) 


NGC4438(S[2) 


NGG4438(F675W) 
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NGC4457(SI2) 




NGC4736(N') 


NGC4736(F555W) 
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NGC5005(N') 


NGC5005(F606W) 


NGC5363(N') 


NGG5363(F300W) 



Figure 2. Gemini MIR and HST optical/UV images of the category I, host-dominated LLAGN with log Lboi/LEdd < —4.6 (no 
HST image is available for NGC 4457). North is up and East left on the images, and a 5" x 5" region is shown. Contours start 
at la and increase in steps of la. All the MIR images were smoothed using a 2-pixel Gaussian. 



NGG4278(N') ' 


• 

NGC4278(F555W) 


• 

NGC4486(N) 


# 

NGG4486(F606W) 


NGC4594(Si2) 


NGC4594(F547M) 





Figure 3. Gemini MIR and HST optical images of the category II, radio-loud LLAGN with log Lboi/LEdd < —4.6. NGC 4374, 
which was not detected in the MIR, also belongs in this category. Contours for NGC 4486 start at lOa and increase in steps of 
lOo", others as in Figure [2] 
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NGC3169(Si2) 


NGC3169(F606W) 


NGC3718(N') 


NGC3718(F606W) 
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NGC3998(SI6) 


NGC3998(F547M) 
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Figure 4. Gemini MIR and HST optical images of the category III, high-Eddington ratio LLAGN (log Lboi/LEdd > —4.6). A 
5" X 5" region is shown for every object except NGC 1097 (24" x 24"). Some of these galaxies were also observed in filters 
other than those shown here, but their morphologies are very similar in all filters. Contours for NGC 1052 and NGC 3718 start 
at 5 a and increase in steps of 5 a, contours for NGC 3998 start ai 7 a and increase in steps of 7 a, those for NGC 3031 and 
NGC 4258 start at 10 a and increase in steps of lOcr, others as in Fig. [2] 
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Figure 5. High-Eddington ratio LLAGN, continued. Contours for NGC5033 start at 5(t and increase in steps of Scr, contours for 
NGC 4579 start at 7cr and increase in steps of 7a, and those for NGC7479 start at 40(J and increase in steps of 40(J. NGC 4261 
is as in Fig. [5] 

4.1.2. The position of the LLAGN on the MIR/X-ray 
,-, relation 
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Figure 6. Detectability of the torus in the LLAGN sample, 
based on the MIR luminosity predicted by the MIR/X-ray re- 
lation (Figure!?]) a-nd assuming a limiting sensitivity of 3 mjy. 
Objects above the solid line are expected to have detectable 
torus emission. Black symbols denote high-Lboi/LEdd ob- 
jects (category III), while radio- loud and host-dominated low- 
Lboi/LBdd nuclei (I, II) are shown by gray and white sym- 
bols, respectively. Galaxies found to have compact MIR mor- 
phologies are given circular symbols, while diffuse galaxies are 
shown with star symbols. NGC 4374, which was not detected 
and whose morphology could therefore not be determined, is 
shown as a triangle. 



can be explained simply by the lower luminosity of the 
AGN relative to the host galaxy. The image morphol- 
ogy alone, then, is not sufficient to demonstrate trends 
of torus properties with Eddington ratio. 



The location of the LLAGN relative to the MIR/X-ray 
relation introduced in the previous section will depend 
on the X-ray data chosen to represent each galaxy. Sev- 
eral X-ray luminosities have been published for some of 
the objects, sometimes with rather different values. This 
was carefully considered during the course of this work, 
and the selection of the intrinsic X-ray luminosities is 



discussed in 

Figure [7] shows that the diffuse, host-dominated 
low-Lboi/LEdd objects, which predominate at the 
low-luminosity end of the sample, tend to have 
higher MIR/lower X-ray luminosities than predicted 
by the MIR/X-ray relation. Four of the six 
host-dominated, low-Lboi/LEdd objects (NGC 4438, 
NGC 445 7, NGC 5005 and NGC 5363) are suspected 
by Gonza lez-Martin et al.l (j2009b[ ) to be Compton-thick 
(Table [5]). This could contribute to an apparent 
MIR excess/X-ray deficit if the intrinsic X-ray lumi- 
nosity is not determined accurately. For the suspected 
Compton-thick nuclei. Table [S] gives intrinsic luminos- 
ity estimates obtained from (a) spectral analysis of 
high-resolution, ~few ke V X-ray data, (b) the ratio of 
Fy(2-10keV)/F([O III]) (IGonzalez-Martm et al.l[2009bl: 
IPanessa et al. 2006"), and (c) the Lboi,AGN — Lfpiy] rela - 
tion (Gouldin g et al.1l2010l:Fereira-Santaella et al.ll2010H . 
The resulting range of estimated intrinsic luminosities for 
the Compton-thick nuclei is indicated in Figure [T] In all 
cases, when the possible large obscuring column is taken 
into account, the AGN could be luminous enough that 
the nuclei are consistent with the MIR/X-ray relation, 
or show a MIR deficit. 

Strong PAH emission (i i4.2p . which lies in the filters 
used to image these galaxies, may enhance the host 
galaxy emission relative to the rest of the sample. It may 
also indicate a contribution from dust heated by circum- 
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nuclear star formation. The PAH emission in the rela- 
tively large-aperture IRS spectra accounts for about 20 - 
40% of the flux in the Michelle/T-ReCS filter bandpasses 
in the host-dominated LLAGN. This is probably an up- 
per limit to the PAH contribution in the ground-based 
images, as the PAH emission is exp ected to be more 



extended than that of the AGN 



e^ 



[2OT0t iHonig et al.|[2010l: lAsmus et all 



Diaz-Santos et ahl 
2011D. The PAH 



emission alone, then, cannot account for the entire MIR 
excess. As noted in §3.11 though, even the combination 
of high-resolution imaging and PSF-scaling photometry 
is unlikely to effectively isolate the weak AGN emission 
in these diffuse objects. Some of the emission must arise 
in the stellar population of the host galaxy; in the pho- 
tospheres of cool stars and in dusty circumstellar shells, 
for example. The contribution of each of these processes 
is difficult to assess reliably with the current data. It 
is possible, however, that once taken into account, the 
AGN MIR emission may be weaker than expected from 
the MIR/X-ray relation, especially if the nuclei are heav- 
ily obscured. 

Finally, we note that emission from X-ray binaries 
may be important in very low-luminosity objects (Lx ^ 
Iq39 gj.g g-i). While this would not explain the excess 
MIR emission, it does further emphasise that the loca- 
tion of very low-luminosity AGN on the MIR/X-ray plot 
may be influenced by a complex mix of factors. 

Of the four radio-loud, low-Lboi/LEdd galaxies, 
NGC 4278 and NGC 4594 have MIR luminosities consis- 
tent with the MIR/X-ray relation, while NGC 4486 has a 
large MIR excess and NGC 4374 (a suspected Compton- 
thick object) has only a MIR upper limit. The images 
of NGC 4278 and NGC 4594 suggest a significant host 
galaxy contribution to the PSF-scaling photometry. Fur- 
thermore, in H4.3.2I we argue that synchrotron radiation 
dominates the MIR emission in these nuclei. Therefore, 
despite the fact that they are ostensibly consistent with 
the MIR/X-ray relation for Seyfert galaxies, it is un- 
likely that NGC 4278 an d NGC 4594 host a Seyfert-like 
torus. We also find (in i j4.3.2p that dust emission con- 
tributes little to the MIR luminosity of NGC 4486. How- 
ever, the MIR/X-ray relation predicts that reprocessed 
IR emission from a torus in NGC 4486 would account for 
< 10% of the observed emission. We cannot estimate the 
amount of nonthermal MIR emission accurately enough 
to rule out torus emission at this level, so taken in iso- 
lation the MIR photometry remains consistent with the 
presence of a torus in this LINER. 

A MIR deficit in the MIR/X-ray plot would be consis- 
tent with suggestions that a Seyfert-like torus does not 
exist in the LLAGN. However, the remaining objects, 
those in the high- luminosity, high-Lboi/LEdd category, lie 
close to the MIR/X-ray relation or exhibit excess MIR 
emission. In most cases the MIR excess is only significant 
at the 1-2(7 level, but it is noteworthy that only two of 
the 13 high-Lboi/LEdd objects lie to the MIR de ficit side 
of the fit. Similar results are reported by ;Asmu s~et al.l 
(|20lH) : while finding that the MIR/X-ray relation is for- 
mally unchanged down to L~ lO^^erg s~^, their LLAGN 
are offset from the relation by about 0.3 dex to the MIR 
e xcess side. 

iHorst et al.l (|2008f ) find that, even in a sample of 
Seyfert galaxies with pointlike nuclei in high-resolution 
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Figure 7. Hard X-ray vs MIR luminosities for the LLAGN and 
for the tiiglier-lumin osity AGN studied by Gandhi ct al. (2009) and 
I Levenson et al.l {2009) . The dashed line shows the fit for the com- 
bined Gandhi and Levenson samples, extrapolated to lower lumi- 
nosities, and the dotted lines mark the +/ — la scatter in the 
Gandhi/Levenson data. Errors on Lx are as stated in the refer- 
ences in Table[5] otherwise assumed to be ±20%. Ellipses show the 
range of available intrinsic X-ray luminosity estimates for possible 
Compton-thick objects, which are identified by encircled markers. 
Other symbols and colors for the LLAGN as in Figure [G] 

MIR imaging, objects with FWHM > 560 x iSub (the 
dust sublimation radius) are systematically offset from 
the MIR/X-ray relation for "well-resolved" AGN. They 
attribute this to contamination from nuclear MIR sources 
other than the torus. All of the LLAGN in this paper 
have FWHM > 560 rsub, so even the high- luminosity, 
high-Lboi/IjEdd with compact nuclei may be subject to 
contamination. The PSF-scaling photometry used in SjS] 
(compared to the aperture photometry of Horst et al.) 
should mitigate this to some extent, but a contribution 
from unresolved, non-torus emission sources - nuclear 
star formation, synchrotron radiation, truncated accre- 
tion disk, etc. - is possible. We return to this issue in 

m 

Finally, NGC 7479 and NGC 3718, two objects with 
compact nuclei and very large MIR excesses, deserve spe- 
cial mention. The very strong silicate absorption feature 
in NGC 7479 fi}42|). unique among this LLAGN sam- 
ple, signifies that this nucleus is deeply e mbedded in a 
conti nuous medium of high optical depth ()Sirockv et al.l 
[2008f) . The nucleus may be Compton thick and the in- 
trinsic AGN luminosity u nderestimated (Tabled i i4.3|) . 
iMufioz Marin et al.l (pOOT f) find the nucleus of NGC 7479 
to be slightly resolved in HST UV imaging, so a nuclear 
star cluster could perhaps also contribute to the dust 
heating and therefore to the MIR luminosity of this ob- 
ject. In the case of NGC 3718, we later present evidence 
that some of the nuclear IR emission is produced by syn- 
chrotron radiation from the jet ( §4.3.2p . 

4.2. 5-35 /im Spectra 

The Spitzer/IRS spectra of the LLAGN are shown in 
Figures |8] - [T3l The area encompassed by the spectra 
(3.6" X 7.2" at 6 /zm) is greater than that shown in 
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Figures H]- [5] It is therefore likely that sources other than 
the AGN also contribute to the IRS spectra, especially 
in the laint, diffuse, host-dominated objects. 

The "bowl-sha ped" spectrum of NG C 4594 (Figure 
[T0|) was noted by IGallimore et all (|2010f). and the other 
radio-loud, log Lboi/LEdd < —4.6 LLAGN have remark- 
ably similar spectra. Silicate emission bands are present 
in at least three of the spectra, including all three type 2 
objects. PAH bands are detected, but are weak in these 
objects. The five host-dominated, low-Lboi/LEdd galax- 
ies with available MIR spectra, on the other hand, all 
show strong PAH emission and (in those objects with 
sufficient wavelength coverage) red continua. 

Spectra are available for nine of the 12 galaxies with 
log Lboi/LEdd > —4.6, and all but one of those galaxies 
(NGC 4261) have compact morphologies in the ground- 
based imaging. Of these nine objects, six show prominent 
silicate emission features. Five of these are type 1.5 - 1.9 
objects, ahhough one (NGC 4261) is a type 2 LINER. 
Strong PAH emission is observed in a further two ob- 
jects (NGC 1097 and NGC 5033, both type 1 objects), 
and the final galaxy in this group, NGC 7479 (SI. 5), has 
very deep silicate absorption bands. NGC 7479 is in fact 
the only object of the 18 in this sample with IRS spectra 
in which the silicate features obviously appear in absorp- 
tion. The continuum properties of these objects are fairly 
diverse. The majority of them have somewhat red con- 
tinua longwards of ~ 25 /xm, but the continuum emission 
of NGC 4261 is essentially flat at A > 20 ^m, and that 
of NGC 1052 is somewhat blue at those wavelengths (in 

The strength and frequency of occurrence of the sil- 
icate emission features in the sample are remarkable. 
For most of the high-Eddington ratio and radio-loud, 
low-Eddington ratio nuclei, we estimate t he strength of 
the 10 /xm silicate feature, 5*10, following ISirockv et al.l 
(|2008D . The resulting values o f Sio are reported in Ta- 
ble [51 IGallimore et al.l (|2010D discuss IRS data reduc- 
tion and analysis at some length, and present Sio for 
a number of ACN based on a full spectral decomposi- 
tion. For the two objects in common with their study, 
we find good agreement in Sio (NGC 4579: 0.40 vs 0.40; 
NGC 4594: 0.31 vs 0.35). Rather than use the simple 
continuum-fitting method to derive for the two log 
Lboi/LEdd > —4.6 galaxies with strong PAH emission 
(N GC 1097 and NGC 50 33), we adopt the values given 
by IGallimore et all ()2010[ ) for those objects. Only three 
of the five host-dominated, low-Lboi/LEdd galaxies have 
full spectral coverage, all have strong PAH bands and the 
Spitzer spectra of these objects probably contain only a 
minor contribution from the AGN. For these reasons, we 
defer measurement of 5*10 in the host-dominated, low- 
Lboi/LEdd galaxies to future work. 

The values of 6*10 in the high- Lboi/LEdd nuclei are 
quite high compared to those commonly observed in type 
1 Seyferts of higher luminosity. Figure [TJ] compares the 
silicate emission strengths in the present LLAGN sample 
and the Seyfert 1-1.5 sample of Thompson et al. (2009; 
also measured using the method of Sirocky et al. 2008). 
The LLAGN sample is neither complete nor unbiased, 
and Figure [T3] is not intended to indicate that the dis- 
tributions of Sio differ systematically between the vari- 
ous AGN types. However, it does demonstrate that the 
10 /xm silicate features found in these particular AGN 



tend to be comparable to the stronger emission features 
in the Thompson et al. Seyfert 1 sample. 

The very pronounced features in NGC 3031 and 
NGC 39 98 have already drawn some attention in the lit- 
erature ([Sturm et al.ll2005l: [Smith et all 120 100. as has a 
further strong feat ure in the type 1 LINER NGC 7213 
([Honig et al.ll2010() . However, strong silicate emission is 
not known to be a general feature of LLAGN. The 10 /im 
emission feature appea rs to be present in t he mean type 
1 LINER spectrum of [Sturm etall ([2006D . but the au- 
thors note that strong silicate emission is not common 
in that sample, which consists of five objects including 
two (NGC 3998 and NGC 4278) in common with this 
paper. As our sample selection was partly based on the 
existence of previous, large-aperture photometry, objects 
with strong silicate emission features could in principle 
have been more likely to be included in this study. As 
the flux increases by less than a factor of two across the 
10 /im feature, though, it is not clear that this would be a 
signiflcant effect. In ij5.3l we suggest that the strong sili- 
cate emission features may represent a particular stage in 
the evolution of the torus. Analysis of the Spitzer spec- 
tra of a larger and/or more targeted sample of LLAGN 
would be helpful in testing that hypothesis. 

4.3. Spectral Energy Distributions 
4.3.1. Data comptlatton 

We combine the new IR photometry with published, 
high-resolution data to produce radio - to - X-ray nuclear 
SEDs for the LLAGN (Table[71 Figures [in]-[2ni)- In terms 
of the nuclear IR emission of LLAGN this SED collection 
is by far the most complete yet available. 

As well as the new MIR data presented here, the SEDs 
include published, ground-based MIR photometry of sev- 
eral sources. These data were obtained using telescopes 
of 5 - 10 m diameter. Ground-based MIR observations 
are close to diffraction-limited under typical atmospheric 
conditions, so these measurements generally have resolu- 
tion <0.5". All the objects for which we include pub- 
lished MIR data are described as having compact nuclei, 
and the photome try is either from small apertures (e.g. 
iHorst et al.ll2008[) or result s from radial profile analysis 
(e.g. [Reunanen et al.|[20Tol ). 

The published 1-5 /xm data are taken from a va- 
riety of sources. The shortest- wavelength points are 
usually from HST FllOW and F160W imaging. The 
flux from the unresolved component is estimated by flt- 
ting a PSF pl us expone ntial, Nuke r or other profile to 
the data fe.g. [Ouillen e t al. 200l|; [Ravindranath et all 
[2OOII : [Alonso-Herrero et al.ii2003i) . thereby decomposing 
the image into a stellar and a nuclear component. Emis- 
sion from stellar photospheres can cer tainly contribute 
in th e NIR in moderate apertures fe.g. [Roche &: AitkenI 
[19851: IWillner et al.l[l985l: Iciemens et al.ll201lD . but the 
unresolved flux given by this procedure should repre- 
sent the AGN alone (unless the nucleus hosts an un- 
resolved stellar cluster). Longwards of 2 /im, little 
HST photometry is available and most of the measure- 
ments we include are from ground-based imaging. Aside 
fro m two galaxies with published adaptive optics imag- 
ing ([Prieto et al.|[2010l: [E^nandez-Ontiveros et al.[[201lL 
Fernandez-Ontiveros priv. communication), the resolu- 
tion of the ground-based data is generally a few tenths of 
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Figure 8. MIR spectra of category I, host-dominated LLAGN with log Lboi/LEdd < —4.6 
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Figure 9. MIR spectra of category I, host-dominated 
LLAGN with log Lboi/LEdd < —4.6, continued. 

an arcsecond, significantly poorer than that of the HST 
imaging. Nonetheless, an unresolved nuclear flux can be 



extracted from the data, or an upper limit given (e.g. 
lAlonso-Herrero et al.|[200l . 

As well as data intended to represent the nuclear emis- 
sion, lower-resolution photometry has also been added 
to illustrate the difference between IR SEDs based on 
large- and small-aper ture measurements. These data 
come from T a ble El lAlonso-Herrero et al.l ()200H) and 
iWillner et al.l (|1985| ) and references therein, and use 
apertures of 2", 3", and 5-7", respectively. 

In constructing the non-IR port i ons of t he SEDs we 
have dr awn o n the work of iHol (|1999D. iNagar et aT 



iMaod ()2007[ ) and, especially. lEracleous et al 



among many others. Additional high-resolution 
data points have also been included if not already present 
in the above compilations, and for objects not covered by 
that work. At all but X-ray frequencies the criterion for 
inclusion is resolution <1". 

For the X-ray measurements, we prefer results based 
on high spatial resolution data. In particular, Chandra 
observations (resolution ^ 1") most effectively spatially 
isolate the nuclear emission and minimise contamina- 
tion from other sources. The effect of pileup in Chan- 
dra observations must be considered for bright sources; 
for instance, the small-aperture luminosity adopted for 
NGC 4486 is based on data taken using short frame times 
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Figure 10. MIR spectra of category II, radio- loud LLAGN with log Lboi/LEdd < —4.6. Dashed lines show continuum fits used 
to estimate silicate feature strength. 

been adjusted to account for this effect. However, Ta- 
ble |6] gives possible intrinsic X-ray luminosities for the 
Compton-thick candidates. 

We have also added some larger-aperture (~7-15") 
submillimeter data to the SEDs. While unlikely to repre- 
sent emission from the nucleus alone, the data provide a 
useful consistency check on estimates of jet emission and 
also give a view of another relatively unexplored spectral 
region. As well as the SED points themselves, Table [7] 
gives the telescope and instrument used for each mea- 
surement, with aperture/resolution information where 
available. The data used in the SEDs were generally 
obtained over a period of several years, and we note that 
variability is known to be significant in so me of these ob- 
jects and in some wave length regimes (e.g. ISchodel et al.l 
.2007tiPian et al.ll20Tol ). 



desig ned to avoid pileup problems (jPerlman fc WilsonI 
[20051) . For sources known to be highly variable, we at- 
tempt to use observations made close in time to the MIR 
observations. Finally, because the derived X-ray proper- 
ties are dependent on the model fitting, where possible 
we choose from results that fully describe the applied 
model and the resulting quality of the fit. 

Some of the objects in this sample are suspected to 
be Compton-thick. In the low-Lboi/LEdd categories, 
IGonzalez-Martm et al.l (j2009a[) use several indicators to 
conclude that NGC 4374, NGC 4438, NGC 4457 and 
NGC 5363 may have Gompton-thick nuclei. High- 
energy Swift Burst Alert telescope observations of 
NGG 7479, a high-Lboi/LEdd object, suggest an in- 
trinsic 2 - 10 keV luminosity several times higher 
than the one adopted in this paper, and modelling of 
this object' s XMM-Newton sp e ctrum implies Nh ^ 
10^-^ cm-^ (iPanessa et all [20061: iCusumano etall [Ml 
iBrightman fc Nandral 1201 lI V Published, extinction- 
corrected 2-10 keV luminosity estimates for these ob- 
jects may therefore underestimate their true intrinsic lu- 
minosities. Because of the large uncertainties inherent 
in attempting to correct for possible Gompton-thickness, 
the SED data in Table [7| and Figures [H] - [101 have not 



4.3.2. Characteristics of the SEDs 

In Figures [TS] - [IHI the SEDs are compared to the 
mean type 1/2 Seyfert SEDs of IPrieto et"all (jlOlO), 
which result from a careful collection of high-resolution 
data spanning a wide frequency range. These aver- 
age SEDs combine data for only 3-4 Seyfert galax- 
ies of each type, so they may not accurately reflect 
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Figure 11. MIR spectra of LLAGN with log Lboi/LEdd > —4.6 (category III). Dashed hnes show conti nuum fits to objects 
without strong PAH emission, used to estimate silicate feature strength. The spectrum of NGC 3998 is from ISturm et al.1 (|2005l ) 
and has had narrow emission features removed. 



the mean properties of the Seyfert population as a 
whole. For example, the MIR/X-ray ratios differ be- 
tween the type 1 and 2 SEDs, contrary to the obser- 
vation that, in larger samples of Seyferts, this ratio in 
the various AGN types is statistically indistinguishable 
iGa ndhi et al. 2009; Lcvcnson et al . 20 09; Honig et al.' 
[201 of) . Nonetheless, the iPrieto et~an 01(X ) Seyfert SEDs 
serve as a valuable comparison for the overall, multi- 
wavelength characteristics of the LLAGN SEDs. A close- 
up view of the 1-20 /im region for each galaxy possessing 
at least two IR data points can be found in Appendix [Xj 
There, the IR SEDs are compa red to the mean 1 - 2 
^m Seyfert SEDs of Ramos Almeida et al.l ()2009l[20Tl.l . 

As expected, the shape of the LLAGN SEDs changes 
markedly in the IR as the photometric aperture decreases 
and the measurements are less contaminated by stellar 
emission. When observed in an aperture of a few arcsec- 
onds, almost all of the SEDs are strongly peaked in the 
NIR. The high-resolution measurements, though, reveal 
a variety of underlying spectral shapes. 

I. Nuclear SEDs of the host-dominated, low- 
Eddington ratio LLAGN 

The SEDs of the host-dominated, low-Eddington ratio 



LLAGN tend to be rather sparsely sampled, with a num- 
ber of non-detections at radio wavelengths. The available 
data show a prominent peak in the MIR but, as discussed 
in fJ3] and the MIR photometry is unlikely to prop- 
erly isolate the nucleus in these objects. The objects with 
high-resolution optical data show rather flat IR-optical 
spectral slopes. In NGC 4736 and NGC 5005, this may 
be due to a contribution from a comp act stellar cluster in 
the NIR/optical (Quillcn et al. 2001: Munoz Marin etHI 
I2007t IGonzalez Delgado et al.il2008ll . 

II. Nuclear SEDs of the radio-loud, low- 
Eddington ratio LLAGN 

The SEDs of these LLAGN are characterised by ex- 
cess radio emission relative to the average Seyfert galaxy 
templates, by the lack of a well-defined MIR peak, and 
by very flat slopes in the MIR - optical region. Type 2 
Seyferts have very steep IR - optical spectra, interpreted 
as extinction by dust in the torus which remains opti- 
cally thick well into the IR. As the hot dust in Seyfert 1 
nuclei is relatively unobscured, their spectra do not show 
the same precipitous drop froin MIR towards NIR wave- 
lengths (lAlon so-Herrero et al.l 120031: IPrieto et all 120101: 
.Ramos Almeida et al.ll20lH) . Despite the fact that 3/4 
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Figure 12. MIR spectra of LLAGN with log Lboi/LEdd > —4.6 (category III), continued. 
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Figure 13. MIR spectra of LLAGN with log Lboi/LEdd > 
—4.6 (category III), continued. 

of the LLAGN in this group are type 2 objects, their 
MIR - NIR/optical SEDs are at least as flat as that of 



the mean type 1 Seyfert (Table [Sj see also Appendix 
for a close-up view of the 1-20 fim region) . 

The four galaxies discussed in this section are among 
the most radio-loud of the present LLAGN sample, so we 
may expect a significant contribution from synchrotron 
radiation to their IR emission. This would have the ef- 
fect of smearing out any underlying, Seyfert-like MIR 
peak in the SEDs and could naturally produce a flat 
MIR - NIR/optical slope. Two of the objects in this 
category (NGC 4374 and NGC 4486) formed part of 
th e sample used to defin e the radio-optical correlation 
of iChiaberge et all ()1999D . who interpret the correlation 
as indicating that the optical luminosity is dominated 
by synchrotron emission. This suggests that (a) a sub- 
stantial IR synchrotron contribution is likely, and (b) 
the core radio and optical SED points may be used to 
assess the magnitude of that contribution. We there- 
fore use parabolic fits to the radio and optical/UV core 
data to determine whether the nuclear IR emission is 
also synchrotron-dominated. This method has been used 
in a number of studies of radio galaxies (iClearv et al.l 
20071: iLeipski et al.ll2009l : Ivan der Wolk et al.ll2010[ ). and 



Landau et al.l ( 198' 



luminous 
such a function 



show that the radio-to-UV SED of 
radio-loud active galaxies can be well fit by 
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Table 5 

IR and X-ray properties of the LLAGN. 



Galaxy Fpsp/F^,,''''^ Sw" "m'*''' log Lmir'' log La^iokcV^ Nh Ref.' CT?g 

(erg s"-"-) (erg s^-"-) (xl0^^cm~^) 



NGC 1052 


0.68 


0.17 


2.4 


41.8 


41.2 


12.90 


1 


N(l) 


NGC 1097 


0.93 


0.21 


1.4 


41.4 


40.6 


0.023 


2 


NGC 3031 


0.97 


0.55 


0.9 


40.5 


40.7 


0.094 


3 


N(2) 


NGC 3166 


0.43 






41.2 


40.7 




4 


NGC 3169 


0.69 






41.2 


41.4 


11.20 


5 




NGC 3718 


0.67 




1.4 


41.1 


40.0 


0.8 


6 




NGC 3998 


0.96 


0.37 




41.3 


41.3 


2.30 


1 


N(l) 


NGC 4111 


0.40 






40.6 


39.6 




7 


N(l) 


NGC 4258 


1.00 


0.30 


2.5 


41.1 


41.0 


6.8 


8 


N(2) 


NGC 4261 


0.60 


0.21 


1.9 


41.2 


41.2 


8.37 


9 


N(l) 


NGC 4278 


0.40 


0.20 


1.3 


39.8 


40.0 


< 0.035 


5 


N(l) 


NGC 4374 




0.19 




< 40.5 


39.6 


0.20 


7 


Y(l) 


NGC 4438 


0.64 






40.9 


39.1 


0.12 


7 


Y(l) 


NGC 4457 


0.50 






40.9 


39.0 


0.1 


7 


Y(l) 


NGC 4486 


0.94 


0.04 


1.5 


41.2 


40.0 


0.14 


10 


N(l) 


NGC 4579 


0.88 


0.40 




41.5 


41.2 


0.93 


1 


N(2) 


NGC 4594 


0.25 


0.31 




39.9 


40.0 


0.19 


1 


N(l) 


NGC 4736 


0.33 






39.8 


38.8 


0.27 


11 


N(l) 


NGC 5005 


0.15 






40.4 


40.1 


3 


12 


?(1) 


NGC 5033 


0.68 


-0.13 


0.8 


41.2 


41.1 


< 0.03 


13 


N(2) 


NGC 5363 


0.27 






40.7 


39.8 


2.66 


1 


Y(l) 


NGC 7479 


0.99 


-2.19 


3.3 


42.7 


41.1 


58 


14 


?(2) 



^ Ratio of MIR PSF-scaling flux to that obtained in a 2" aperture. See ^3.11 

^ Based on photometry in the Si-2 (8.8 ^m) filter where available, otherwise N' (11.2 fim), Si-5 (NGC 1097; 11.7 
fim) or N (NGC 4486; 10.8 fim). See Table|2]and El 

Strength of the ~ 10 fira silicate feature ( i|4.2ll . Positive values of 5io indicate emission. Typical errors on Sio 
are estimated to be about 0.05. 

NIR - M IR s pectral index, oc u~°'if-, for objects possessing at least one MIR and one H band photometric 
point. See ^431 

" Absorption-corrected 2-10 keV luminosity. See ^4.31 for details of X-ray data selection 

* Reference for X -ray luminosity and col umn density. 1: [Gonzalez-M artm et al.l ll2009bl V 2: INemmen et al.l H2006I V 
3: 'Swartz et al. ("2003). 4: derived from Fabbiano et al. (1992) using standard assumptions to convert from count 
rate to flux. 5: Terashima & Wilson (2003). 6: Satyapal et al. (2005). 7: Flohic et al. (2006). 8: Young & Wilsoij 
(|2004|) . 9: Zezas et al. (2005). 10: Pcrlman & Wilson (2005). 11: .Eracleous et al.l (120021) . 12: iGuainazzi et all 
(I2TO ). 13: Cappi et al. (2006). 14: Pancssa et al. (2006). 

s Combination of X-ray and optical indicators suggests galaxy may be Compton-thick. The hx and Nh values 
in colum ns 5 and 6 have not bee n modified to account for possible Compton-thickness. 1: IGonzalez-Martfn et al.l 
(I2009al) . 2: IPanessa et al.l 1 I2OOI) . 



Table 6 

Intrinsic X-ray luminosity estimates for suspected 
Compton-thick LLAGN. 



Galaxy logLx(spec.)- logLx([0 III])'' logLx([0 IV])'^ 





(erg s 1) 


(erg s 1) 


(ergs 1) 


NGC 4374 


39.6 


41.3 


39.3 


NGC 4438 


39.1 


40.8 


40.7 


NGC 4457 


39.0 


40.6 


40.4 


NGC 5005 


40.1 


41.6 


40.6 


NGC 5363 


39.8 


41.6 




NGC 7479 


41.1 


42.9 


41.6 



As given in Table |5] from spectral fitting of high-resolution X- 
ray data. Corrected for absorption, but will nonetheless under- 
estimate the intrinsic X-ray luminosity of a truly Compton-thick 
AGN. 

^ " Corrected" lum i nositi es from lGonzalez-Marti'n et al.l (|2009bl ) 
and IPanessa et al.l (|200^ ), based on Fx (2 - 10 keV)/F([OIII]) 

ratios 

Based on the L]-,oi a.gn ~ L[oiv] relation of 'Goulding et al.l 
(I2OIOI), using published [O IV] fluxes (Dudik et al. 2009; 
IPereira-Santaella et al.ll2010l : iPanuzzo et al.ii2011 ) and assuming 

Lboi/Lx = 16 dHaEosg, E} 



As shown in Figure [TTl the fits indicate that in ah four 
cases the small-scale IR emission is dominated by the 
jet. This was already proposed for NGC 4486 (M87) by 
iPerlman et al.l (1200 If). It is also consistent with the re- 
sults of lvan der Wolk et all ([2010') , who find that most of 
the FR-I and low-excitation FR-II objects in their sample 
have IVIIR emission consistent with a synchrotron spec- 
trum. In the case of NGC 4486 we reiterate that, al- 
though the vast majority of the IMIR emission is likely to 
originate in the jet, the expected luminosity of the torus 
in this object, predicted by the MIR/X-ray relation, is 
only a small fraction of the observed MIR luminosity 
(Figure ^ . The uncertainty in fitting the IR jet contri- 
bution is very likely large enough to allow the existence 
of such a faint torus, so the SED fitting alone does not 
rule out the presence of a Seyfert-like torus in NGC 4486 
and similar radio galaxies. The lack of a steep IVIIR - 
NIR/optical spectral slope is consistent with an unob- 
scured nucleus, but could also be explained by, for in- 
stance, significant synchrotron emission on scales larger 
than a ~pc-scale torus. 

On the other hand, the observed IVIIR emission of 
NGC 4594 is already consistent with the standard 
IVIIR/X-ray relation, and subtracting the estimated syn- 
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Table 7 

Spectral Energy Distribution Data 



Galaxy 


V 


vLv 


Aperture / resolution^ 


Facility 


Reference 




(Hz) 


(erg s-i ) 


arcsoc 






NGC 1052 


5.00E09 


5.43E+39 


0.002 


VLBA 


Maoz (2007); Kadler et al. (2004) 


NGC 1052 


8.40E09 


9.04E+39 


0.002 


VLBA 


Maoz (2007); Kadler et al. (2004) 


NGC 1052 


1.50E10 


1.41E+40 


0.0005 


VLBA 


Kovalev et al. (2005) 


NGC 1052 


2.20E10 


1.49E+40 


0.0005 


VLBA 


Kadler et al. (2004) 


NGC 1052 


4.30E10 


1.29E+40 


0.0003 


VLBA 


Maoz (2007); Kadler et al. (2004) 


NGC 1052 


8.60E10 


1.07E+40 


0.001 


VLBI 


Lee et al. (2008) 



^ Aperture diameter, resolution or beam size, as 
To be published as an online-only table 
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Figure 14. Strength of the 10 /im silicate emission 5*10 
in the high-Lboi/LEdd (category III, gray) and radio-loud 
low-Lboi/IjEdd (category II, cross-hatched ) LLAGN com- 
pared with Sio in the Seyfert 1 sample of [Thompson et al.l 
|2009). For clarity, NGC 7479 from the LLAGN sample and 
UGC 5101 from the Seyfert 1 sample, which have deep ab- 
sorption features (S'lo = -2.19 and -1.52 respectively), are not 
shown. 

chrotron component would create a significant MIR 
deficit relative to the expected torus emission in this ob- 
ject. NGC 4594 is therefore likely to be a genuine un- 
obscured type 2 AGN. In all of these objects, however, 
even if a weakly-emitting torus does exist, its contribu- 
tion to the total nuclear MIR luminosity will be negli- 
gible for the purpose of even fairly detailed modeling of 
the broadband SED. 

III. Nuclear SEDs of the high-Eddington ratio 
LLAGN 

The SEDs of the high-Eddington ratio LLAGN are 
quite diverse, more so than those of the radio-loud, low- 
Eddington ratio objects. If they have one distinguishing 
feature, it is that many of them are more radio-loud than 
the average Seyfert galaxy. This is expected from previ- 
ous work showing that ra dio loudness increases with de- 
creasing Eddington ratio (jHcill2003 iTerashima fc WilsonI 
|2003; Sikoraet al. 2007|). 

Of the better-sampled SEDs, some - those of 
NGG 1052, NGC 4258, NGC 4261, NGC 4579 and 
NGC 7479 - show a well-defined, Seyfert-like peak at 
MIR wavelengths. The MIR - NIR/optical spectral 
slopes in these objects are within the range brack- 



3ed in or inferred from the original paper. 

eted by the iPrieto et all ()2010D Seyfert 1 and 2 tem- 
plates. Amongst these objects, nuclei of the same 
AGN type do not necessarily have the same spectral 
slope; for instance, the type 1.9 Seyferts NGC 4258 and 
NGC 4579. However, if intermediate-type LLAGN have 
similar nuclei to more luminous intermediate-type ob- 
jects, then a variety o f spect ral shapes is to be expected; 
lAlonso-Herrero et al.l ()2003l ) find a wide range of spec- 
tral indices in Seyferts of types 1.8 and 1.9 in the CfA 
sample. Judged solely by their SEDs, there is little indi- 
cation that these objects are anything but "scaled-down" 
Seyferts, complete with a MIR "dust bump", albeit with 
extra radio emission. The high-Eddington ratio objects 
in this sample are almost exclusively type 1-1.9 nu- 
clei, so we cannot search for the systematic differences in 
SEDs, spectra etc. that are characteristic of type 1 and 
2 Seyferts of higher luminosity. 

A few of the galaxies - NGC 1097, NGC 3031, 
NGC 3718, NGG 3998 and NGC 5033 - have rather 
flat MIR - NIR/optical SEDs compared to the Seyfert 
templates, and no well-defined MIR peak. These were 
also characteristics of the radio-loud, low-Eddington ra- 
tio objects, whose IR emission in i i4.3.2l we found to 
be dominated by the jets. In NGC 1097, the "excess" 
NIR/optical emission is probably related to the young, 
mass ive nuclear star duster (Storchi-Bcrgmann et al] 
120051 : [Mason et al.i 120071 ). The strong PAH bands and 
extended 3 /xm emission in NGC 5033 may indicate the 

S iresence of nuclear star form ation in that object as well 
Alonso-Herrero et al.l I2003D . For the remaining three 
galaxies, we use the fits described in i)4.3.2l to investigate 
the possible IR synchrotron emission. The fits suggest 
that, although synchrotron emission is unlikely to domi- 
nate the MIR in the same way as for the lower-luminosity, 
radio-loud, low-Eddington ratio LLAGN, by contributing 
in the NIR it could flatten the spectral slopes as observed. 
This fitting procedure assumes that the measured opti- 
cal luminosities represent optical jet emission, however, 
which may not be a valid assumption for these objects. 

As an alternative means of investigating whether the 
flat MIR - optical SED slopes are indicative of signifi- 
cant IR synchrotron emission, we searched for a correla- 
tion between the IR spectral index (Table [5]) and radio 
loudness (log vL^ih GHz) / Lx, © for those galax- 
ies in i i4.3.2l and in this section which have at least one 
high-resolution NIR data point. We found no relation 
between these quantities. However, there are several pos- 
sible ways in which such a correlation could be masked. 
As well as emission from compact clusters such as in 
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Figure 15. SEDs of host-dominated LLAGN with log Lboi/LEdd < —4.6 (category I). Open circles denote data with resolu- 
tion/aperture > 1". The X-ray luminosity predicted by the MIR/X-ray re lation ( t |4.1.2|l for each galaxy is shown as an open 
square. Solid and dashed lines mark the mean type 1 and 2 Seyfert SEDs of IPrieto ct al. (2010), respectively, normalised in the 
MIR. 

NGC 1097, extinction in some objects may mask an in- teristics of these LLAGN. 
trinsically flat slope. A case in point is NGC 4261, one 
of the more radio-loud objects in the sample whose nu- 
cleus is obscured by Nh ~ 10 cm" "2 (,Zezas et al..:2005l) 
and which is suspected to be heavil y extinguished in the 
optical/UV ()Eracleous et al.ll2010bf ). The lack of a corre- 
lation could simply reflect that a diversity of phenomena 
governs the properties of the LLAGN SEDs, and IR ob- 
servations of a larger sample of objects would be useful 
in this respect. 

To summarise, many of the high-Eddington ratio 
LLAGN have SEDs that broadly resemble those of con- 
ventional Seyferts, although often with enhanced radio 
emission. In some cases, though, the MIR - NIR/optical 
SED slopes are flatter than expected for Seyferts, which 
hints at possible stellar and/or jet contributions in the 
IR. Other sources, such as a truncated accretion disk, 
may also influence the IR properties of the SEDs. De- 
tailed modeling of the SEDs with jet, RIAF, thin disk 
and torus components is beyond the scope of this paper, 
but in ij5.3l we suggest a scenario which is consistent with 
the morphological, spectral and broadband SED charac- 
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Figure 17. SEDs of radio-loud LLAGN with log Lboi/LEdd < —4.6 (category II). Dot-dashed lines show parabolic fits to the 
nuclear radio/optical/UV data ( ii4.3.2|) . Other lines and symbols as in Figure [T^ 
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Figure 18. SEDs of high-Eddington ratio LLAGN (log Lboi/LEdd > —4.6, category III). Symbols and lines as in Figure [TTl 
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Figure 19. SEDs of high-Eddington ratio LLAGN (category III), continued. 
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Figure 20. SEDs of high-Eddington ratio LLAGN (category 
III), continued. 
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5. DISCUSSION 

We have compiled subarcsecond-resolution, 1-20 /im 
imaging of 22 LLAGN, including new IR observations of 
20 objects. The data have been used to investigate the 
morphology of the objects and their place on the stan- 
dard MIR/X-ray plot, and to fill in a region of the SED 
hitherto lacking high-resolution information. In addi- 
tion, we have presented Spitzer IR spectroscopy of 18/22 
objects. 

As a pilot study intending to provide an overview of 
the IR properties of LLAGN, the sample is somewhat 
heterogeneous and may be biased towards objects bright 
in the MIR. Rather than discuss statistical trends within 
the sample, then, in this section we highlight some in- 
teresting properties of the galaxies, loosely divided into 
three groups with certain common elements. 

5.1. /. The host-dominated, low-Eddington ratio nuclei 

At the low-luminosity end of the sample are nuclei with 
very low Eddington ratios (log Lboi/LEdd < —4.6) that 
are faint compared to the surrounding host galaxy emis- 
sion. If a dusty torus is present in these objects, it is 
very likely too faint to contribute significantly to even the 
high-resolution, PSF-scaling photometric measurements 
(Figure [5]) . 

All of the host-dominated, low-Eddington ratio galax- 
ies exhibit strong PAH emission in their central re- 
gions. Although commonly used to identify starburst 
galaxies, PAH bands are also observed in many other 
environments suc h as the diffuse I SM of dusty el- 
liptical galaxies (iKaneda et al.l I2008D . Three of the 
five LLAGN in the host-dominated, low-Eddington ra- 
tio category (NGC 4438, NGC 4457 and NGC 5005) 
are included in the ste llar population analysis of 
iCid Fernandes et"an (|2004[ ) and lGonzalez Delgado et all 
()2004[) . Intermediate-age stars (10^ — 10^ yr) contribute 
~30 - 60% of the optical light in these and many of the 
other LLAGN in those studies, but the contribution of 
stars aged < 10'' yr is ^ 5% or less. The PAH emission 
in these objects may be related to the intermediate-age 
stellar population; based on analysis o f the band ratios in 
a small sample of early-type galaxies, iVega et all (|20100 
conclude that the PAH molecules are supplied by mass- 
losing carbon stars formed within the last few Gyr. Fur- 
ther detailed and quantitative studies of the PAH bands 
in LLAGN would be of value in determining whether the 
features are related to the known intermediate-age pop- 
ulation or whether they reveal active star formation that 
is difficult to detect in conventional optical studies. 

5.2. //. The radio-loud, low-Eddington ratio nuclei 

Also at the low-Eddington ratio, low-luminosity end of 
the sample, we highlight a set of galaxies with very strong 
radio emission. The nuclear SEDs of these objects do 
not show the prominent mid-IR peak observed in "con- 
ventional" Seyfert galaxies, and their MIIR-NIR/optical 
spectral slopes are flatter than that of the average type 
1 Seyfert. We find that the nuclear IR emission is 
dominated by synchrotron radiation from the jet, con- 
sistent with results for various radio galax y samples 
(|Leipski et al.|[2009Hvan der Wolk et al.ll2010D . 

High-resolution SEDs can show whether the IR emis- 
sion is energetically dominated by nonthermal processes. 



However, additional information is needed to conclu- 
sively rule out the presence of a Seyfert-like torus. The 
well-known MIR/X-ray relation can be used to estimate 
the MIR emission that would be expected from the torus, 
and in at least one galaxy (NGC 4594) the observed 
MIR luminosity is already no more than would be ex- 
pected from a torus, before accounting for the likely 
dominant MIR synchrotron component. This galaxy is 
therefore probably an example of a genuine, "bare" type 
2 AGN with no obscuring torus (also suggested by the 
low X- ray column, 2 x 10^^ cm^^; Gonzalez- Martin et al.l 
l2009bD . The same may be true for NGG 4486. How- 
ever, the emission from a torus heated by the weak AGN 
in NGC 4486 would account for only -10% of the ob- 
served luminosity, well within the error of the SED fit- 
ting used to investigate the synchrotron component. In 
terms of the IR emission there is therefore room to "hide" 
a torus in NGC 4486. Additional evidence, such as 
the minimal X-ray absorption towards this typ e 2 nu- 
cleus ()Di Matteo et al.l[200l iPerlman fc Wilsoni[2005l ) is 
needed to support any assertion that NGC 4486 is an 
unobscured type 2 AGN. 

The IR spectra of the galaxies in this category, 3/4 
of which are type 2 objects, all appear to show silicate 
emission. The PAH and [Ne II] emission in NGC 4278 
mean that the presence of the silicate features in this 
object is ambiguous, and the 11.3 /im PAH band in 
NGC 4374 could contribute significantly to the measured 
value of Sio in that nucleus. However, the PAH fea- 
tures in NGC 4486 and NGC 4594 are we ak and the sili- 
cate emission in NGC 4 486 is well-known (jPerlman et al.l 
120071 iBuson et al.ll2009t ). Silicate dust features are an un- 
ambiguous sign that dust is present. In higher-luminosity 
AGN, silicate emission features are characteristic of type 

1 objects and are an expected res ult of a direct v iew of 
hot dust in a roughly face-on torus (jNenkova et al.| [2008') . 
Silicate emission is also known in a few typ e 2 AGN and 
can be reproduced by cluinpy torus models (iMason et al.l 
120091: iNikutta et alT 120091 : lAlonso-Herrero et al.l 120 lit ). 
This explanation seems unlikely for the radio-loud, low- 
Eddington ratio nuclei for two reason s. First, silicat e 
emission is rare in type 2 Seyferts (e.g. iShi et al1l2006D . 
so it would be surprising to encounter it in all three type 

2 objects in this sample. Second, as noted above, some 
of these galaxies probably do not host a torus. 

If not arising in a standard, Seyfert-like torus, the sili- 
cate emission bands could arise in diffuse, optically thin 
dust perhaps associated with the remains of a dissipating 
torus. We develop this line of reasoning in H5.3[ Alter- 
natively, the features may be produced in circumstellar 
dust shells. The weak, extended silicate emission in sev- 
eral early-type ga laxies has been show n to originate in 
mass- losing stars (jBressan et al.|[2006[ ). and this process 
is also likely to exp lain the silicate features in NGC 4486 
(|Buson et al.l[20Ci9l ). 

5.3. ///. The high- Eddington ratio nuclei 

The nuclei with log Lboi/LEdd > —4.6, which lie at 
the high luminosity end of the sample, tend to have 
prominent nuclear point sources in the MIR. They of- 
ten have strong silicate emission features compared to 
those typically observed in Seyfert 1 nuclei, and their 
MIR luminosities are consistent with or slightly in excess 
of those predicted by the standard MIR/X-ray relation. 
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The nuclear broadband SEDs of these objects are rather 
mixed. Some are essentially indistinguishable from those 
of "conventional" Seyferts, many have excess radio emis- 
sion compared to higher-luminosity Seyferts, and a few 
have unusually flat MIR - NIR/optical slopes. 

Considered in isolation, these characteristics offer no 
compelling evidence that the torus is absent in these 
LLAGN or that it differs from that observed in higher- 
luminosity Seyferts. However, we do find some indica- 
tions that the IR emission may not arise in the standard, 
Seyfert-like torus of the AGN unified model. In Figure 
\n\ we plot the strength of the 10 /im silicate feature 
ag ainst HI c olumn density, for both the Seyfert galaxies 
of iShi et al.l ([2006) and the LLAGN with bio measure- 
ments (fg^S Table E]). In the Shi et al. Seyferts, Sio 
and Nh are loosely correlated. The scatter presumably 
reflects the fact that the observed 5io is likely a com- 
plicated funct ion of the precise arran gement of clouds in 
the torus C Honig fc Kishimotoll201Cl[) . and may also in- 
dicate a contribution from absorption in the host galaxy 
(iRoche et al.ll2007HDeo et al.l2009HAlonso-Herrero et al.1 
120 111) . Maiiv of the high-Eddington ratio LLAGN lie on 
the upper envelope of the Seyfert points, with relatively 
high values of Sio per unit N//. One possible explana- 
tion is that the dust-to-gas ratio is lower in these par- 
ticular objects than in most Seyfert galaxies. This may 
be expected if the torus is an optically thick region of 
an accretion dis k wind that become s less powerful at low 
accretion rates (Elitzm & Ho 2009) , when the amount of 
material reaching the dust sublimation radius and able to 
form dust grains is reduced. In this case we may expect 
either a torus with fewer clouds and a higher probabil- 
ity of observing a hot, directly-illuminated cloud face, or 
simply optically thin dust emission. Both of these con- 
figurations would cause relatively strong silicate emission 
features, consistent with our finding that Sio in many of 
these high-Eddington ratio LLAGN is strong compared 
to that in typical ty pe 1 Seyferts (Figure [HI) . 

As mentioned in a contribution to the observed 
silicate features from circumstellar shells in the host 
galaxy is also possible. As the luminosity of the central 
engine diminishes, the relative strength of features aris- 
ing in the surrounding stellar population will increase. 
This could also account for changes in 5io/Nh among 
LLAGN. 

The high-Eddington ratio LLAGN have at least as 
much MIR emission as predicted by the MIR/X-ray re- 
lation for Seyferts and quasars. If the torus in these 
LLAGN does indeed contain less dust, the IR continuum 
emission must be produced by some other mechanism. 
A contribution from synchrotron radiation is likely, and 
may explain the unusually flat MIR -NIR/optical slopes 
observed in some of the SEDs. Another possibility is that 
some of the IR emission comes from a truncated accretion 
disk. In the models of iNemmen et all ()2011[ ). the disk 
emission peaks between 1-10 /xm and can account for 
essentially all the luminosity at these wavelengths. Emis- 
sion associated with nuclear st ar clusters, such as that 
known to exist in NGC 1097 (iStorchi-Bergmann et aLl 
l2005[) may also play a role. However. I Asmus et al.l (|2011f l 
find that star formation contributes <30% of the small- 
scale 12 jim flux in most LLAGN of comparable lumi- 
nosity to those studied in this paper. Detailed modelling 
of the LLAGN, taking advantage of the new constraints 



from the well-sampled IR SEDs presented here, will pro- 
vide more insight into the processes responsible for the 
IR emission of LLAGN and the role of the torus in these 
objects. 



• High-Eddington ratio LLAGN 
eRadio-ioud, iow-Eddington ratio LLAGN 
□ Seyfert (Shi et ai. 2006) 
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Figure 21. Strength of the 10 /im sihcate feature vs X-ray 
column density for the high-Eddington ratio (category III) 
and radio- loud, low-Eddington ratio (category II) LLAGN. 
Square symbols indicate the Seyfert galaxies of IShi et aLl 
((200Q1 . The dashed line shows Shi et al.'s fits to their Seyferts, 
the dotted line their fit to their whole sample (mostly com- 
prised of Seyferts and various classes of quasar) . Positive val- 
ues of Si o in dicate emission, negative values absorption. As 
in Figure lT4l NGC 7479, which has a deep absorption feature 
(Sio = —2.19), is omitted for clarity. 



6. CONCLUSIONS 

The new, high-resolution IR imaging and SEDs pre- 
sented in this paper represent the first step towards es- 
tablishing the nuclear IR properties of a significant num- 
ber of LLAGN. At the median distance of the galax- 
ies studied here, 16.8 Mpc, the AGN is well isolated in 
objects with log L2_io keV > 40.5 erg s~^ (log Lboi > 
41.8 erg s^^). Deeper observations would allow a cleaner 
detection of the pointlike central engine in less luminous 
or more distant objects, but the cost in observing time 
would be large with current facilities. 

The Spitzer spectra of the less radio-loud of the lowest- 
luminosity LLAGN show strong PAH features. Detailed 
analysis of the band ratios etc. may demonstrate a con- 
nection with very young stars or active star formation, 
rather than the well-known intermediate-age stellar pop- 
ulation or the diffuse ISM. If so, combining MIR and 
optical spectroscopy of a larger sample could provide a 
more complete view of the star formation histories of 
LLAGN. 

Some of the lower-luminosity, low-Eddington ratio 
galaxies are very strong radio emitters. The high- 
resolution SEDs of these objects lack a well-defined MIR 
peak and are very flat in the MIR - NIR/optical range. 
The broadband SEDs allow us to conclude that their 
IR emission is dominated by synchrotron radiation from 
jets. Some of these objects, such as NGC 4594 (the Som- 
brero galaxy) are probably unobscured type 2 objects, 
genuinely lacking a broad line region. In other objects. 
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such as NGC 4486 (M87), the torus emission predicted 
by the standard AGN MIR/X-ray relation is so weak that 
a torus cannot be ruled out based on SED information 
alone; supporting information is also needed. 

At relatively high Eddington ratios and luminosities, 
the galaxies have prominent, compact MIR nuclei. Their 
SEDs are quite diverse, ranging from objects indistin- 
guishable from typical type 2 Seyferts to nuclei with 
strong radio emission and flat MIR - NIR/optical slopes. 
The strong silicate emission features present in many of 
these objects demonstrate that dust is present in their 
nuclear regions. However, we tentatively suggest that 
those strong features and the values of Sio/Njy in these 
nuclei may indicate optically thin dust and low dust-to- 
gas ratios, consistent with models predicting the disap- 
pearance of the torus in LLAGN. If so, much of the IR 
continuum emission must be caused by processes other 
than emission by dust in the torus. Possible sources in- 
clude synchrotron radiation, a truncated accretion disk 
and nuclear star formation. Observationally, IR po- 
larimetry and monitoring campaigns may constrain the 
role of each of these processes. A targeted study of the 
silicate features in a well-defined LLAGN sample should 
also help clarify the properties of the torus in these ob- 
jects. Finally, we expect that detailed modelling of the 
data in terms of accretion disk, RIAF, jets and torus will 



be valuable in elucidating the nature of the IR emission 
of LLAGN. 
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APPENDIX 

APPENDIX 1: THE 1 - 20 AtM REGION OF THE SEDS 

The high-resolution IR data are the novel aspect of this paper, so i n this Appendix we present figures focusing on 
the 1-20 fim part of the SEDs. Rather than the mean Seyfert SEDs of Pricto ct al. (2010) used in H4.31 in this s ection 
the LLAGN SEDs are compared to the mean SEDs of RaffiosAlmcida c t al, (2009, 2011). The Prieto etall (|20100 
templates have the advantage of wider wavelength coverage, but those of lRamos Almeida et all (|2009l 1201 Ih " include 
more objects. We therefore judge that the Ramos Almeida data are more suitable for a detailed comparison of the 1 
- 20 ^m region. Both sets of SEDs are shown in Figure [22l While the type 1 SEDs agree fairly well (and would be 
indistinguishable on the scale used in Figures [T5] - BU]) . the type 2 SED of (f'rieto et al.l ()201 [)l) is significantly less steep 
than that of Ramos Almeida ct al. (2009). The difference in the type 2 IR SEDs does not affect the discussion in ij4.31 
which covers a broad wavelength range and highlights differences between some of the LLAGN an d the average type 
1 Seyfert. IR spectral indices for these LLAGN are given in Table [3 For comparison. iRamos Almeida et aLr(|2011l ) 
give mean values of 1.6 ± 0.2 for type 1 Seyferts and 3.6 ± 0.8 for type 2 Seyferts. 
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Figure 22. 1 - 20 fim SEDs of the radio-loud, low-Eddington ra tio (category II) LLAGN hav ing at least two IR photometric 
points. The solid black line shows the mean type 1 Seyfert SED of lRamos Almeida et al.l (|2011l '). the dashed black line the mean 
type 2 Seyfert SED of Ramos Almeida et al. (200^). T he LLAGN are d enoted by the dotted line and filled circle symbols. For 
NGC 4278 only, the average Seyfert 1 and 2 SEDs of Pricto et al.l (|2010l ) are plotted in gray. The template SEDs are normalised 
at 8.7 fim. The observed SEDs are normalised either directly t o the Si2 8.7 /im) photom etric point, if available, otherwise 
using the filter closest in wavelength and assuming qmir = 2.0 (E amos Almeida et al.ll201ll V Errors on the observed data are 
not shown as error estimates are not available for many of the photometric points from the literature. Where available, errors 
are typically comparable to the size of the plot symbols. 
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Figure 23. 1 - 20 /im SEDs of the high-Eddington ratio (category III) LLAGN having at least two IR photometric points. 
Lines and symbols as in Figure 
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Figure 24. 1 - 20 /^m SEDs of the high-Eddington ratio (category III) LLAGN having at least two IR photometric points, 
continued. 
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